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We study the formation of local magnetic moments in quantum point contacts. Using a Hubbard- 
like model to describe point contacts formed in a two dimensional system, we calculate the magnetic 
moment using the unrestricted Hartree approximation. We analyze different type of potentials 
to define the point contact, for a simple square potential we calculate a phase diagram in the 
parameter space (Coulomb repulsion - gate voltage). We also present an analytical calculation of 
the susceptibility to give explicit conditions for the occurrence of a local moment, we present a 
simple scaling argument to analyze how the stability of the magnetic moment depends on the point 
contact dimensions. 

PACS numbers: 73.23.-b,75.75.+a 



I. INTRODUCTION 

The problem of electric charge transport through quan- 
tum point contacts (QPC) has been intensively studied 
both from the theoretical and the experimental points 
of viewA The observation of conductance quantization 
in a diversity of point contacts is now well established 
and, in general terms, understood. The discovery of ex- 
tra structure, that looks like conductance plateaus at 
0.7(2e 2 /fo) in GaAs devices, however, still remains as 
an open questionAi4Afiiii2i£ The recent accumulation of 
experimental evidence as well as the theoretical analy- 
sis of this structure suggest that it is due to magnetic 
fluctuations. 10 ! 11 ! 12 

The experimental data around the anomalous 0.7 
plateau observed in some of these QPC have been in- 
terpreted as due to a Kondo effect This interpreta- 
tion is based on the observation of a zero bias anomaly 
in the conductance, its temperature and magnetic field 
dependence as well as a single energy scale ksTx as- 
sociated with it. The Kondo effect is due to the mag- 
netic screening of a localized spin, a phenomenon that 
occurs for magnetic impurities diluted in metalsA^ and 
also for quantum dots in mesoscopic circuits While all 
the 0.7 anomaly features observed in QPC are consistent 
with the occurrence of Kondo effect, it is not clear how a 
magnetic moment could develop in these systems. Point 
contacts in GaAs-AlGaAs heterostructures are built as 
narrow constrictions in a two dimensional electron gas 
formed at the interface between GaAs and AlGaAs. The 
constriction is built using patterned surface depletion 
gates so the shape and size of the QPC can be controlled. 

It has been shown that a two dimensional electron gas 
with a constriction under some special conditions, that 
can be achieved by tuning a gate voltage, can develop a 
magnetic moment. This has been done by using spin de- 
pendent density functional theory 15 on the one hand and 
unrestricted Hartree-Fock solutions of effective Hubbard 
models on the other M As already pointed out, one should 



be careful in the physical interpretation of the frozen 
spin solution obtained with these methods. As in the 
old impurity problem of magnetic moment formation^ 
the mean field solution can not be correct since it breaks 
the local symmetry. In the exact solution of the problem 
spin fluctuations should recover local rotational invari- 
ance. It is very difficult to incorporate spin fluctuations 
on top of the mean field solution to fully recover the rota- 
tional invar iance. However the mean field solution gives 
a good indication of the region in parameter space where 
we should expect magnetic fluctuations to play a central 
role in the low energy physics™ What is still missing in 
the problem of the magnetic nature of QPC is a detailed 
analysis of the condition for the occurrence of a magnetic 
moment including its geometry and size dependence. In 
this work we use the Hartree criteria to determine the re- 
gion of the parameter space where the contact develops a 
magnetic moment. In next section we present the model 
and its mean-field version. In section Hi Al the numerical 
solution is used to determine the region of stability of a 
localized magnetic moment in the point contact. In part 
III Bl the analytical expressions for the mean field mag- 
netic instability are presented, and in part III Cl the limit 
of narrow resonances is used to predict size scaling. The 
last section includes summary and discussion. 



II. THE MAGNETIC INSTABILITY AT THE 
POINT CONTACT 

A natural description of the two dimensional (2D) elec- 
tron gas in GaAs- AlGaAs heterostructures is an effective 
mass theory in which the kinetic energy is given by the 
one of a 2D Fermi gas of particles with an effective mass 
m* and a characteristic particle density n ~ 
For this system the Fermi wave vector kp is of the order 
of 10 6 /cra. The potential created by the applying gate 
voltages is described by a function V(r) where r is the 2D 
coordinate. Finally, to describe the magnetic properties 
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of the system, the electron - electron interaction has to 
be included explicitly. In order to solve the Schrodinger 
equation in the Hartree approximation, we discretize the 
space to end up with an effective tight binding model with 
hopping matrix element t = h 2 /2m* a 2 where a is the dis- 
cretization parameter, the lattice parameter of the tight 
binding model. The potential energy V(r) is included 
as an on site energy and the Coulomb repulsion as an 
on site Hubbard term. The discretization parameter a 
should be taken such that the Fermi wave vector kp is 
much smaller than any reciprocal lattice vector, typically 
kp <C 27r/a. 

Our starting point is a two dimensional extended Hub- 
bard Hamiltonian 17 with a constriction and a gate poten- 
tial, H = Ho + Hi n t with 

#0 = ^2^4^ ~ t J2 C t C 3v (!) 
i,(T (ij),cr 

and 

Him = U[J2 C i^ C t C il + \ J2 Vij4* C iv C ^jvi (2) 

Here cf a creates an electron at site i with spin a, the 
diagonal energy e% is a function of the coordinate of site 
i and defines the point contact with the gate potential. 
The last term in Hq the the site independent nearest 
neighbor hopping. The parameter U is the local Coulomb 
repulsion and the dimentionless parameter rjij < 1 gives 
the spacial dependence of a screened interaction. 

In the Hartree approximation the potential energy is 
renormalized by the Coulomb repulsion 

H M F = ^SiaC^Civ - t ^ C t* C j* ~ K, (3) 
i,a (ij),cr 

with s i(T = Si + U[(ni-o) + Ej,^ 1 ~ s ij)Vij ( n j>')]> 
ni,a = c t(j c iv ls the number operator and (..) indicates 
the expectation value and K is a constant. In what fol- 
lows we use the notation 

n i = ( n i^) + ( n hi) 

m % = (K T > - (m d ))/2 (4) 
so that Sicy = Si — crU rrii with = Si + U[rti/2 + ^^(1 — 

The shape of the QPC as well as the gate voltage is de- 
fined through the potential energy Si. We study systems 
with the geometry of Fig. ^ consisting of a two dimen- 
sional strip infinitely long and with a width of N sites. 
The QPC is defined as a constriction at the center of the 
strip. In the rest of the work, we assume a square lattice 
and weak Coulomb repulsion U, much smaller than the 
Stoner critical value. We also take the average charge 
per site no < 1 to avoid any Fermi surface nesting ef- 
fect. Then far from the QPC the charge is uniform and 
the magnetization is zero so that the potential energy is 
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FIG. 1: a) Schematic picture of the model for the QPC (dark 
region) between the source and drain leads, the potential at 
the constriction is V g . b) Profile of the potential along the 
QPC. 



£icr — £7no [1/2 + $^(1 — Sij)rjij]. To do the calculation, 
we artificially divide the sample in two regions: a cen- 
tral region (between vertical dashed lines in Fig. [Q that 
includes the QPC and the uniform regions that include 
the source and drain far away from the point contact. In 
the numerical calculation we evaluate the self-consistent 
solution within the unrestricted Hartree scheme for the 
central region coupled with the lateral regions with uni- 
form charge and zero magnetization. This procedure is 
acceptable if the charge and magnetization profiles are 
continuous and smooth at the boundary between regions. 
The same scheme is used to analyze the spin susceptibil- 
ity. 

A. Numerical Solution of the Unrestricted Hartree 
Equations 

In this section we present the numerical solution of 
the unrestricted Hartree approximation. The charge rii 
and magnetization rrii at each site of the central region 
containing the QPC are evaluated with the solution of 
the self-consistent equations: 

M = -\f F duIm[G^(w)] (5) 

here the retarded Green function Gf a {u) is the diagonal 
element of 

G R = [lj + z0+ - H° MF - Z R (u + z0+)] _1 (6) 

where H^ F is the Hartree Hamiltonian of the central 
part and S jR (o; + zO + ) is the self energy due to the lateral 
regions with rn a = 0. 

We studied different QPC shapes. In the present ap- 
proximation, the long range part of the electron-electron 
interaction renormalizes the on site energy and redefines 
the shape of the QPC. From hereon we take rjij = 0. We 
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FIG. 2: (a) Magnetic moment (thin line) and total charge 
(thick line) at the constriction of a 3-site-wide and 9-site-long 
QPC. The local interaction is U — 2.0 and the Fermi energy is 
€f — —1, with all parameters in units of t = 1. The source and 
drain slabs have a width of 20 sites, b) Local density of states 
at the Fermi level, averaged over the sites of the constriction, 
calculated in the restricted Hartree approximation. 




FIG. 3: Spatial distribution of the a) magnetization, c) charge 
for V g = 0.66, and other parameters as in Fig. |2 b) and d) 
Same as a) and c) respectively with V g = 2.27. 



first present results for a model square potential shown 
in Fig. [I] and denned as: Si = 0, V g and oo for i in the 
source and drain, in the neck of the QPC and at the sides 
of it respectively. We considered a variety of QPC with 
No = 3 sites in the lateral direction and N\ = 9 sites 
in the longitudinal direction. We define the total charge 
and the magnetization of the contact as Q — rii and 

M = J2 i mi where the sum is over all the sites of the 
QPC. The self-consistent results are shown in Fig. Ofor a 
QPC of width No and length N±. For large gate voltages 
the total charge Q in the point contact is exponentially 
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FIG. 4: Phase diagram U—V g for a rectangular 3x9 QPC, dark 
regions indicate higher values. Left: magnetization at the 
QPC, Right: negative derivative of the charge with respect 

tO Vg. 



small. As V g decreases Q increases and for some values of 
the gate potential there is an abrupt increase in the total 
charge. The steps obtained around this point correspond 
to approximately one electron been transferred from the 
source and drain to the point contact. This behavior in Q 
is characteristic of Coulomb blockade. Between the two 
first steps in Q there is a spin 1/2 localized at the point 
contact as indicated by the magnetization curve versus 
V g shown in Fig. Efa). The local density of states at the 
QPC shown in Fig. EJb) presents a series of resonances. 
The resonances are associated with longitudinal modes 
in the QPC. The occurrence of a local magnetization at 
the QPC coincides with a narrow resonance crossing the 
Fermi energy. The first resonance has no nodes in the 
QPC as shown in the magnetization profile of Fig. Efa). 
As V g decreases other resonances cross the Fermi energy 
and new magnetic solutions are obtained. When the first 
resonance of the second channel cross the Fermi energy, 
again a spin 1/2 is localized in the QPC. The wave func- 
tions of the second channel have a node at the center 
of the QPC in the transverse direction and generate the 
magnetization profile shown in Fig. Efb). For the value 
of U considered in Fig. [^magnetic moments are obtained 
only for the first resonance of each transverse channels. 
However, the other wider resonances may also produce 
local magnetic moments for larger values of U. In Fig. 
0|a) we present the phase diagram in the parameter space 
[U— V g ] for the 3x9 QPC presented up to now. The 
dark regions correspond to a stable magnetic moment at 
the QPC. Fig. |3fb) illustrates the behavior of the charge 
Q; to analyze the behavior of Q it is convenient to plot 
its derivative —dQ/dV g . For small values of U we ob- 
tain non-magnetic solutions for all values of V g . However 
as V g increases and the QPC resonances cross the Fermi 
level, Q increases and there is a maximum in —dQ/dV g . 
For the parameters of the figure, magnetic solutions are 
obtained for U/t > 1 and for values of the gate poten- 
tial V g that make the resonances to coincide with the 
Fermi level. The occurrence of a magnetic solution is ac- 
companied by a splitting of the (—dQ/dV g ) maximum, 
a characteristic of the Coulomb blockade regime. Fig. 
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FIG. 5: Phase diagram U-V g for a rectangular 5x6 QPC, dark 
regions indicate higher values. Left: magnetization at the 
QPC, Right: negative derivative of the charge with respect 

tO Vn. 



shows the phase diagram and the behavior of Q for a 
shorter an wider point contact. In this case larger values 
of U are required to produce magnetic moments. In next 
section we interpret these results in terms of the suscep- 
tibility evaluated in a simple approximation and make a 
scaling analysis to describe how the magnetic instability 
depends on the QPC size. The square potential optimizes 
resonances at the QPC and consequently favors the for- 
mation of local moments. We end this section showing 
some results obtained with a smoother and more realistic 
potential defined as: 



Si ~- 

{1 



VgfiVi ~ 

■ 10 [/Or; 



■oi/2)[l 
+ ao/2) - 



- f(Vi 



■ai/2)] x 
+ a /2)]} ,(7) 



where = (xi,yi) is the coordinate of site z, f(x) = 
(1 + e x / A ) -1 with A a characteristic screening length, ao 
and a\ are the QPC width and length respectively. To 
describe a wedge-like point contact we take a width that 
varies linearly with yf. ao(yi) = ao + 5ao x (yi + a\/2)/a\. 
The potential profiles for these contacts are shown in Fig. 
|H1 In the same figure the average density of states at the 
point contact is also shown. For a wire-like point contact 
resonant states are obtained, the energies of these reso- 
nances are at the bottom of each channel. The density 
of states has much less structure in the case of wedge- 
like contacts. This shows that wire-like contacts are good 
candidates to develop magnetic fluctuations while wedge- 
like structures are not. 



B. Spin Susceptibility 

In the presence of an external magnetic field, the local 
magnetization is given by 



(8) 



where hj is the magnetic field at site j and Xi,j is the non- 
local susceptibility. In our system with no translational 
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FIG. 6: a) Average local density of states in the QPC, at the 
Fermi level for the rectangular smooth potential of c); [7 = 0, 
V g — 0.5, £f — — 1, ao = 3, a\ = 15. b) Same as a) for the 
wedge shaped potential of d) and V g = 1.25, 5ao = 4.5. 



symmetry the non-local susceptibility depends on the two 
coordinates and rj. In the Hartree approximation, the 
energy shift of the one particle levels is used to define 
an effective magnetic field given by hj = h ext + Urrij 
where we assume a uniform external field h ext . The 
magnetization is then given by 
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U 7 \Xi,j m 3 +Xihext 



(9) 



with Xi — ^2jXi,j- As we show below, x% is J us ^ ^ ne 
local density of states at the Fermi energy and for U = 
the Pauli susceptibility depends on the coordinate as 
the local density of states. For h ext = a non-trivial 
solution of equation (0 gives the onset of a spontaneous 
magnetization. We define the susceptibility matrix 
with matrix elements Xi, j •> an d a magnetization vector M 
as a column vector with components m^. Then equation 
(0 for h ex t — has the form of an eigenvalue problem 



(10) 



as U — > there is no non-trivial solution of this equa- 
tion. The onset of a spontaneous magnetization is given 
by the largest eigenvalue of x being equal to 1/U , the 
corresponding eigenvector gives the magnetic profile of 
the instability. 

Since in the paramagnetic state the system has spin 
rotational invar iance, we calculate the transverse suscep- 
tibility. The linear response of the system to a magnetic 
field along the x direction is: 



Xi. 



(ii) 
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where indicates the retarded Green function. 

In terms of the self-consistent one-particle eigenstates of 
the Hamiltonian Hhf with wave functions and en- 

ergies e v , the susceptibility is given by 



VjJL 



(12) 



where f(e) is the Fermi function. 

Due to the orthogonality of the one-particle wave func- 
tions, at zero temperature, we have: 



Xi 



E^ = -Ei^)i 2 ^ ( 13 ) 

3 v V 

\ipv(i)\ 2 5(e u - e F ) = pi(e F ) 



where we have taken the ratio between the Fermi function 
difference and the energy deference as the Fermi func- 
tion derivative when e v — > £ M , p%{e) is the local density 
of states and e F the Fermi energy. Then, inserting the 
above expression in equation JjJ), for U = we obtain 

For U =^ in a system with no translational invariance 
we have to calculate the non-local susceptibility, a ma- 
trix of dimension equal to the number of sites of the sam- 
ple. To simplify the problem we assume that U is much 
smaller than the Stoner value that generates a global in- 
stability, consequently a non trivial solution of equation 
should have the magnetization concentrated in the 
region of the point contact. We assume that far from the 
contact rrii = and look for solutions with mi ^ only if 
i belongs to a small region R that includes the contact. 
Then we have to calculate a reduced matrix susceptibility 
Xi,j with ij in R. 

Since even in the paramagnetic state, and due to the 
charge redistribution close to the point contact, the one- 
particle states of Hh F depend on U in a non-trivial way, 
the largest eigenvalue n of the susceptibility depends on 
U and the instability condition is a self-consistent equa- 
tion: 

k{U) = 1/U 

and its solutions have to be obtained numerically. 



C. The Resonant State Approximation 

Here we present some analytical results based on the 
fact that, for high gate voltages, the local density of states 
at the point contact present well defined resonances. By 
varying the gate voltage the position of the resonances 
can be tuned to coincide with the Fermi level. When this 
occurs, for small quantum point contacts where the quan- 
tization effects are important, the transport and mag- 
netic properties of the contact are dominated by a single 



resonant state and in what follows we consider this situ- 
ation. This is a valid approximation as long as the width 
of the resonance remains much smaller that the separa- 
tion between resonances. The non-local susceptibility is 
now given by 

xij = -E^(^(*K(jK(j) /(£ ; ) : f M (14) 

here the decorated sum indicates summation over all 
states belonging to a single resonance. For i in the point 
contact, the corresponding wavefunctions can be written 

as 



MOI 2 = l«™«l 2 



{e v - A) 2 + 7 2 



(15) 



where a m (i) is the wave function of the QPC, 7 and A 
are the width and the energy of the resonant state m and 
p v is the density of states. Since we can work with real 
wavefunctions, they can be taken as the square root of 
the above expression and the susceptibility can be put as 

Xi,j = Wm(i)?Wm(j)\ 2 Xres (16) 

with the susceptibility of a resonant state given by 

(7A) 2 



Xr 



& ^[(e,-A)2 + 7 2][( £p -A) 2 +7 2 



(17) 



The susceptibility matrix has the form x — XresA 
where the matrix A can be put as 



i« m (i)i 2 

l« m (2)| 2 



[|a m (l)| 2 ,|a m (2)| 2 ..|a m (iV)| 2 ] (18) 



|a m (iV)| 2 

and the instability condition becomes 

aUXres 1 



(19) 



where a is the largest eigenvalue of the matrix A. Prom 
the form of A it is clear that the vector 



M = 



K(i)l 2 

K(2)| 2 
L |a m (iV)| s 



(20) 



is an eigenvector with eigenvalue a = |«m(^)| 4 and 
that all other eigenvalues are zero. The condition for the 
formation of a magnetic moment at the point contact is 



effXres 



(21) 



where U e ff — Si l a m(^)| 4 ^ is the effective Coulomb re- 
pulsion for two electrons at the QPC state m with spacial 
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TABLE I: Critical values of U for the appearance of magnetic 
solutions. 



N x Ni 


U c /t 


u;/t 


3 x 11 


0.55 


0.57 


3x9 


0.85 


0.84 


4x7 


1.50 


1.76 


5x6 


2.25 


2.49 


5x5 


3.18 


3.96 



wavefunction a m (i). If the resonance is centered at the 
Fermi level, this condition is simply 



eff 



7T7/2 



(22) 



For small U and an arbitrary potential form of the 
point contact defined by the potential Vfa), the Hartree 
solution can be used to estimate a m (i) and 7. 

Now we compare this criterion with the full unre- 
stricted Hartree calculation for the lowest energy reso- 
nance. For the square potential of Fig. [I] the U = 
eigenfunctions ao(i) are: 



a (i) 



■ sin( 



(iVo + 1) y a(N + l) 



■ sin( 



+ v a(A^i + l) 



Xi) x 



Vi) 



(23) 



here a is the lattice parameter, clNq and aNi are the QPC 
width and length respectively. As this wavefunction is 
hybridized with the right and left reservoirs it acquires a 
width 7 = 27rpVeff wnere P is the density of states of the 
reservoir and the effective hybridization is 

" 1 



t 



iEedge 



(24) 



the sum is over all the sites of the QPC that are at one 
edge (right or left), hybridized with the reservoir. With 
this estimation and the condition of equation ((221) we 
obtain the critical value U* of the Coulomb repulsion 
for the occurrence of a magnetic solution shown in table 
1. The comparison with the values obtained using the 
fully unrestricted Hartree approximation is very good in 
particular for long and narrow point contacts. 

For the second longitudinal resonance of the first chan- 
nel in the case of the 3x9 geometry, we find the same 
U e ff as for the first resonance while the width of the 
resonance becomes 7' ~ 3.67 [as observed in Fig. Efb)], 
therefore giving U'Jt ~ 3.1 in good agreement with the 
corresponding phase diagram shown in Fig. 21 

Finally we can make a approximate scaling analysis for 
large systems: the effective repulsion in a single resonance 
Ueff = Ei \am(i)\ 4 U ~ U/(N x Ni) and the effective 
width of the resonance is proportional to the square of 
the hybridization of equation (|24j) . 7 ~ pt 2 /Nf. These 
estimations lead to a critical value of U that at reso- 
nance scales with the size of the point contact as U c ~ 



pt 2 No/N±. Numerical estimations of the mean field crit- 
ical value for large systems are in agreement with this 
scaling. 



III. CONCLUSIONS 

We have presented results for the formation of lo- 
cal magnetic moments in point contacts. We used a 
Hubbard-like model to describe point contacts formed 
in a two dimensional system. The contact is defined 
in terms of a potential V(r) that can be varied with 
a single parameter V g representing a gate voltage. We 
calculate the magnetic moment using the unrestricted 
Hartree approximation. For a square potential, the sys- 
tem shows a marked tendency to form a localized mo- 
ment at the point contact each time a new channel is 
tuned to the Fermi energy. In this conditions the crit- 
ical value of the local repulsion U is almost an order 
of magnitude smaller than the Stoner critical value for 
an homogeneous system. For the parameters of figure 21 
the critical value for the first resonances is U ~ t. Us- 
ing the effective mass and electron density characteristic 
of GaAs-AlGaAs hetherostructures, we can take a lattice 
parameter a = lOnra. With these numbers, the results of 
the figure correspond to a point contact of 30nm x 90nm 
and the critical value of U gives an effective interaction 
Ueff ~0.5meV. 

In long contacts defined with a square potential, the 
second longitudinal resonance may also generate a lo- 
cal moment for moderate values of U. This is a conse- 
quence of the square potential that optimizes resonances 
each time the Fermi wavelength is commensurate with 
the contact length. For more realistic potentials only the 
first longitudinal resonance of each channel may gener- 
ate a local moment. Moreover, for wedge-like contacts 
we found no evidence of moment formation in the Hub- 
bard type models. The numerical results are interpreted 
in terms of a simple one-resonance approximation. We 
also present a simple a scaling argument to interpret the 
general dependence of the magnetic instability with the 
point contact dimensions. 

We end by stressing that the Hartree calculation, that 
breaks the spin symmetry, only gives a criterion that 
allows to identify the region of parameter space were 
the low temperature physics may be dominated by mag- 
netic fluctuations. In this particular regions a Kondo-like 
model may be used to describe the spin fluctuations-Sii^. 
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